home *** CD-ROM | disk | FTP | other *** search
/ Collection of Tools & Utilities / Collection of Tools and Utilities.iso / ada / gnat1792.zip / gnat179b / t-adainc / s-xp_bml.adb < prev    next >
Text File  |  1994-05-19  |  3KB  |  55 lines
  1. ------------------------------------------------------------------------------
  2. --                                                                          --
  3. --                         GNAT RUNTIME COMPONENTS                          --
  4. --                                                                          --
  5. --                        S Y S T E M . X P _ B M L                         --
  6. --                                                                          --
  7. --                                 B o d y                                  --
  8. --                                                                          --
  9. --                            $Revision: 1.1 $                              --
  10. --                                                                          --
  11. --           Copyright (c) 1992,1993,1994 NYU, All Rights Reserved          --
  12. --                                                                          --
  13. -- GNAT is free software;  you can  redistribute it  and/or modify it under --
  14. -- terms of the  GNU General Public License as published  by the Free Soft- --
  15. -- ware  Foundation;  either version 2,  or (at your option) any later ver- --
  16. -- sion.  GNAT is distributed in the hope that it will be useful, but WITH- --
  17. -- OUT ANY WARRANTY;  without even the  implied warranty of MERCHANTABILITY --
  18. -- or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License --
  19. -- for  more details.  You should have  received  a copy of the GNU General --
  20. -- Public License  distributed with GNAT;  see file COPYING.  If not, write --
  21. -- to the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. --
  22. --                                                                          --
  23. ------------------------------------------------------------------------------
  24.  
  25. function System.Xp_BML
  26.   (Left  : Long_Long_Unsigned;
  27.    Right : Natural)
  28.    return  Long_Long_Unsigned
  29. is
  30.    Result : Long_Long_Unsigned := 1;
  31.    Factor : Long_Long_Unsigned := Left;
  32.    Exp    : Natural := Right;
  33.  
  34. begin
  35.    --  We use the standard logarithmic approach, Exp gets shifted right
  36.    --  testing successive low order bits and Factor is the value of the
  37.    --  base raised to the next power of 2.
  38.  
  39.    --  Note: it is not worth special casing the cases of base values -1,0,+1
  40.    --  since the expander does this when the base is a literal, and other
  41.    --  cases will be extremely rare.
  42.  
  43.    while Exp /= 0 loop
  44.       if Exp rem 2 /= 0 then
  45.          Result := Result * Factor;
  46.       end if;
  47.  
  48.       Factor := Factor * Factor;
  49.       Exp := Exp / 2;
  50.    end loop;
  51.  
  52.    return Result;
  53.  
  54. end System.Xp_BML;
  55.